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We discuss a relativistic chiral theory of nuclear matter with a and lu exchange using a 
formulation of the a model in which all the chiral constraints are automatically fulfilled. We 
establish a relation between the nuclear response to the scalar field and the QCD one which 
includes the nucleonic parts. It allows a comparison between nuclear and QCD information. 
Going beyond the mean field approach we introduce the effects of the pion loops supple- 
mented by the short-range interaction. The corresponding Landau-Migdal parameters are 
taken from spin-isospin physics results. The parameters linked to the scalar meson exchange 
are extracted from lattice QCD results. These inputs lead to a reasonable description of the 
saturation properties, illustrating the link between QCD and nuclear physics. We also derive 
from the corresponding equation of state the density dependence of the quark condensate 
and of the QCD susceptibilities. 
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§1. Introduction 

o. 

t^. , The subject of this work is the interplay between nuclear physics and QCD and 

examples of what each field can bring to the other. Its motivations have been twofold. 
The first aim is to study QCD related quantities such as the quark condensate in the 
nuclear medium or the QCD scalar susceptibility which is its derivative with respect 
to the quark mass, in a way which is fully consistent with the saturation properties 
of nuclear matter. The second one is to build a relativistic theory of nuclear matter 
which satisfies all the chiral constraints. It is possible to find a common frame to 
reach this goal in effective theories. They are built to mimick QCD at low energy, 
^ ■ allowing the study of QCD related quantities. Moreover a linear realization of chiral 

symmetry involves a scalar isoscalar field which can generate the nuclear attraction, 
i.e., the scalar field of the aoj model of Walecka-Serot. 1 - 1 Our final conclusion will be 
that there exists a direct and model independent link between QCD quantities and 
the parameters which govern the attraction in the atu model. In the linear sigma 
model the explicit symmetry beaking part of the Lagrangian is C x sb = ca where 
c = f n m 2 and a is the scalar field, chiral partner of the pion. This quantity plays 
the role of the symmetry breaking Lagrangian of QCD : C^gg = -2m q qq. We 
have then the equivalence : 

qq(x) _ a(x) 



(qq)vac U 
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Thus it is the expectation value of the scalar field which controls the evalution of 
the condensate : 

(qq(p)) = ( a (p)) / 1-2 n 
(qq)vac U 

and the sigma propagator (taken at zero momentum) which controls the susceptibil- 
ity : 

Xs = J d 4 x (qq(x) qq(0)) = Mp£ J d 4 x(a(x) a(0)} = Mf^O). (1-3) 

For the nuclear binding problem we use, as in ref., 2 ) a non linear representation, 
keeping a chiral singlet scalar field which is the radius of the chiral circle, S. In the 
vacuum S = f w , and we denote s the fluctuation: S = f n + s. In the usual non linear 
sigma model this fluctuation is ignored. We have shown previously, (ref. 2 )) that the 
mean value of s can be identified with the scalar mean field of the au model. In this 
way, with the introduction of the non-linear representation, all chiral constraints are 
fulfilled. This would not be the case if the identification would be made instead with 
the non chiral invariant sigma field. The two fields s and a are related through the 
transformation from polar to cartesian coordinates : 

a = (f n + s) cos f(J-^J ~U + s-^- (1-4) 

in which terms of order s cj) 2 have been ignored. The effective nucleon mass is influ- 
enced by the mean field s : 

M* N (s) =M N (l + -f) ~ M N - 4r PS- (1-5) 

Thus the nucleon mass and the condensate evolutions are linked but not proportional. 
The term which involves the scalar density of nuclear pions 4> 2 is absent in the mass 
evolution. Similarly the propagator of the a field, which gives the scalar susceptibility 
and the one of the s field are linked but not identical, with (ref. 3 )) : 

D* = D S + ^G (1-6) 

where G is the two-pion propagator. At low momenta the s field is not coupled to two 
pions in contradistinction with the a one which is strongly coupled (said differently 
the s field, which is at the origin of the nuclear binding, just decouples from low 
energy pions the dynamics of which is described by chiral perturbation theory). 

In this framework a joint study of the nuclear binding and of the QCD quantities 
is natural and the compatibility between the two is insured. The question is wether 
something more can be learned from this common study. In a first step we ignore 
the pion loops. In this case the condensate evolution is simply given by : 

A (qq(p)) = J_ , 1-7 x 

(qq)vac U 
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Fig. 1. Sigma self-energy in the nuclear medium. 



where the symbol A for the condensate stands for the modification with respect to 
the vacuum value. Similarly the nuclar susceptibility, vacuum value subtracted, is : 



Xs 



2 

vac 



f 2 

J 71 



(D*-D° 



(1-f 



Here D* is the full in-medium sigma propagator while D® is the vacuum one. Both 
are taken at zero momentum and they are related by : 



(D*-D Q 



*2 



to* 2 



n s (o) 



*2 ' 



(1-9) 



In this expression, g$ is the scalar coupling constant and 17$ is the full RPA particle- 
hole polarization propagator in the scalar-isoscalar channel. The mass to* is the 
effective sigma mass which differs from the free one mostly by the effect of the 
tadpole term : 



,*2 



1 + 



3s 

7, 



+ 




3ffs 



■ps 



(MO) 



This represents a large decrease of the mass ~ 30% at po. The tadpole contribution 
and the effect of the p.h. polarization propagator are illustrated in fig. [H The first 
part corresponds to the a self-energy from the following o~N (non Born) amplitude, 
T a N = —3g s /fn, while the second part arises from the in medium modified (Pauli 
blocked) o~N Born amplitude. Inserting the expression (|1-9|) of the sigma propagator 
into that, eq. (jl-8p . of the susceptibility we obtain to lowest order in density : 



Xs 



{qq) 



2 

vac 



3gg ps 



+ g s n s (o) 



Ml) 



In this expression the first term linear in the density embeds the scalar susceptibility 
of the individual nucleons while the second term reflects the effect of the nuclear 
excitations which affect the QCD susceptibility through the coupling of the quark 
density fluctuations to the nucleonic ones. 4 ) All in all the parenthesis on the r.h.s. 
of eq. (|l-lip is the total nuclear response to the nuclear scalar field, 1Z . The 



4 



M. Ericson and G. Chanfray 



proportionnality factor between 1Z A and Xs can be expressed in terms of the scalar 
quark number of the nucleon arising from the a field : 

qs = = JJSh^l f d 3 x(N\a(x)\N) = -M^l MS (142) 
2 m, Jit J fa mi 

such that : 

X A s = 2^Tl A . (1-13) 

9s 

It is remarkable that the full response is reflected in the QCD one (and vice versa), 
including the nucleon part. 6 ) 

Coming back to eq. (jl-lip the term linear in the density provides the QCD 
scalar susceptibility of the nucleon from the scalar field : 

(*?)* = "2^^- (M4) 

Jit '"■a 

This component, which to our knowledge has not been signalled before, is an example 
of a new information that the nuclear problem can bring. The question is wether 
this component has a reality. In the model certainly since it can be also be obtained 
as the derivative of Q s s with respect to the quark mass. And in reality? There 
are indications in favor of a negative component (beside the pionic one discussed in 
ref. 5 )) in the lattice results on the evolution of the nucleon mass with the quark 
mass, (equivalently the pion squared mass). They are available only above ~ 
400 MeV and in order to extract the physical nucleon mass an extrapolation has 
to be performed. Thomas et al. 7 ) have separated out the pion which introduces a 
non analytical behavior in m q . For this they evaluate the pionic part of the self- 
energy, U n , introducing different shapes of the form factor at the ttNN vertex with 
an adjustable cut-off parameter. They expand the rest of the nucleon mass in powers 
of m 2 as follows : 

M N (ml) = a + a 2 ml + a 4 m^ + E n (m n ). (1-15) 

The parameters dj show little sensitivity to the shape of the form factor, with values 
d2 — 1.5 GeV^ 1 while 04 ~ —0.5GeV~ s , (see ref. 7 )). From this we can infer the 
non-pionic pieces of the sigma commutator using the Feynman-Hellman theorem : 

a non- P ion = = ^2 + 2ai m* ~29 MeV . (1-16) 

It is dominated by the a 2 term. Its value indicates the existence of a large component 
in o"tv beside the pion one, that it is natural to attribute to the scalar field. The next 
derivative provides the non pionic susceptibility : 

(ria) 2 8 / non-pion\ , v 2 
non— pion \HHIvac I N 1 \ L l L l/vac a /l T7\ 

Xns ~ ^Tdmi\~^T~)~^r A - { } 



It has a negative sign, as expected from the scalar contribution. If the signs are 
right, are the numerical values also compatible with the sigma model? If we identify 
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,i , non—mon j non—pion -,i ,1 j 1 1 a i. 

the two quantities, a N and Xns with the sigma model values, we hrst 

have : 

non—pion 2 s s 2 9S r -. -. Q \ 

a N " ~ a 2 m n = a N = f^m^^- (1-18) 

which leads to m a = 800 MeV (for 55 = Mj^/f^). On the other hand the identifi- 
cation of the two expressions 11-141 and 11-171 for x^~ p%on and the elimination of m a 
using [T48l leads to the relation : 

a 4 ~ — ^-a 2 . ~ -3.4 GeV~ z (149) 
2M 2 v ; 

while the value found in the expansion is only —0.5 GeV~ ?J . The model as such fails 
to pass the QCD test. In fact this is to be expected and even gratifying because 
it also fails the nuclear physics test. Indeed the effect of the tadpole term on the 
a propagator is too large. The softening of the sigma mass makes nuclear matter 
collapse and prevents saturation. 8 ) We have seen that the QCD susceptibility and 
the aN amplitude are related. Both values are too large in magnitude in the model 
which is incomplete and has to be improved. Indeed an important effect is missing, 
namely the scalar response of the nucleon, k^vs, to the scalar nuclear field, which 
is the basis of the quark-meson coupling model, (QMC), introduced in ref. 9 ) The 
crucial point is that its origin lies in confinement in such a way that its sign is 
positive, i.e,, it opposes an increase of the scalar field. In QMC the response is 
calculated in a pure bag model. The quantity kns is then related to the QCD scalar 
susceptibility of the bag through a relation similar to ours (jl-13p but in which only 
bag quantities enter. The positive sign follows as the confined quarks become less 
relativistic with increasing quark mass and their scalar number increases. In order 
to incorporate all aspects, it will be interesting to explore the full nucleonic response 
to the scalar field in a model of the nucleon, such as the one of Shen and Toki, 10 ) 
which incorporates both aspects, i.e., where the nucleon mass originates from both 
the coupling to the condensate and the confinement. 

It is possible 11 ^ to improve our previous model described above with the phe- 
nomenological introduction of the scalar nucleon response, kns- It modifies the 
nucleon mass evolution as follows, 

M* N (s) = M N (l + jj + i k ns s 2 . (1-20) 

The sigma mass is also affected : 

mf~ml -(^- K NS )PS, (1-21) 

Jit 

and the nucleonic QCD susceptibility as well, in such a way that the relation between 
04 and 02 becomes : 

fl4 = "2M (3 " 2C) - (1 ' 22) 
where C is the dimensionless parameter C = (f%/2M) k^s- Numerically 04 = 
— 0.5GeV~ 5 gives C = +1.25. A large cancellation of the tadpole effect is indeed 
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required by the lattice expansion, total cancellation occuring for C = 1.5. Our ap- 
proach has then consisted in going from QCD to nuclear physics and use the lattice 
expansion to fix the parameters of the ouj model. 12 ) The question is wether this 
procedure makes sense and wether it leads to a possible description of the nuclear 
binding. We have explored this, first in the mean field approach, but also with the 
introduction of the pion loops, 12 ) as is described below. In the latter case the con- 
densate acquires a new component linked to the nuclear pion scalar density, (0 2 ). 
However the nucleon mass is blind to this component. 2 ) This is also true for the 
sigma mass and the eq. (|l-20p and (jl-2ip remain valid. Pion loops nevertheless 
affect the energy. They do not enter at the mean field level but contribute through 
the Fock term and through the correlation energy schematically depicted in fig. [2j 
It includes iterated pion exchange and also the part of the NN potential from the 
two-pion exchange with A excitation. In the correlation term we have introduced be- 
side pion exchange the short-range components via the Landau-Migdal parameters, 
g' NN , g' NA , g' AA . We have taken their values from a systematic survey of the data on 
spin-isospin physics by Ichimura et al. 13 ^ For the parameters g' they indicate large 
deviations from universality with : g' NN = 0.7, g' NA = 0.3, g' AA = 0.5. Beside pion 
exchange we have also introduced the transverse channel, dominated by p exchange 
together with the short-range component. The other parameters of the model are 
chosen as follows. The form factor at the irNN vertex is taken as a dipole with a 
cutoff parameter A = 7 , which is guided by the following considerations : it leads 
to a pionic contribution to <jn of 22 MeV . Adding the lattice value of the non-pionic 
part, a non - pion = 29 MeV, the sum takes the value a N = hi MeV, in the well ac- 
cepted range. The scalar coupling constant is the one of the model gs = Mf f n = 10 
and for the sigma mass we have followed the lattice indications, allowing a small 
readjustment around the lattice value, which is m a = 800 MeV. We have found 
a better fit with m a = 850 MeV which corresponds to a non - pion = 26 MeV. The 
omega mass is known and the loNN coupling constant totally free. For the nucleon 
scalar response we have followed the indications of the lattice data but not strictly in 
view of the uncertainties attached to the higher derivatives. The value which fits the 
saturation properties is found to be C ~ 1, not far from the lattice value C = 1.25. 
With these inputs we have obtained a satisfactory description of the nuclear binding. 
The binding energy per particle is shown in fig. [3l with its different components. We 
like to comment on the correlation energy which has a value of — 17 '.4 MeV . The 
longitudinal channel is strongly suppressed by the short-range component, (with less 
cancellation in its NA component). This leaves the transverse part as the domi- 
nant contribution (— 9.5 MeV versus —7.9 MeV for the longitudinal one). From 
the equation of state and using the Feynman-Hellmann theorem we have deduced 
the quark condensate. The next derivative with respect to m q provides the scalar 
susceptibility. The two QCD susceptibilities, scalar and pseudoscalar are depicted 
in fig [U . The second one follows the evolution of the condensate 4 ) and it is remark- 
ably linear in the density in spite of the interaction. The scalar one shows a strong 
increase with respect to the free value, surpassing even the pseudoscalar one around 
2po- This behavior is largely due to the mixing with the nuclear excitations. 
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Fig. 2. Schematic representation of the longitudinal spin-isospin contribution to the correlation 
energy. 
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Fig. 3. Binding energy of nuclear matter with = 8.0, m a = 850 MeV and C = 0.985. The 
full line corresponds to the full result, the dotted line represents the binding energy without 
the Fock and correlation energies and the dot-dahed line corresponds to the contribution of the 
Fock terms. The decreasing dotted line (always negative) represents the correlation energy. 



In summary we have studied a chiral relativistic theory of nuclear matter based 
on the sigma model and its interplay with QCD. We have shown that the soften- 
ing of the sigma mass arising from the tadpole term, (from 3cr interaction) has a 
counterpart in QCD in the form of a negative contribution to the nucleonic QCD 
scalar susceptibility from the scalar field. We have investigated the presence of this 
component in the lattice expansion of the nucleon mass with respect to the quark 
mass. There is indeed an indication in favor of a negative component but its magni- 
tude is much too small. This is in fact totally consistent with nuclear physics. The 
only effect of the tadpole with the softening of the a mass prevents saturation. In 
both cases a canceling effect must occur. For this we have introduced, as in QMC, 
the scalar response of the nucleon which is a reflect of the QCD scalar susceptibility 
of the nucleon which is due to confinement. Our approach has then be to utilize 
the QCD information on the nucleon mass evolution with the quark mass to fix or 
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1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

0.8 1 1.2 1.4 1.6 1.8 2 2.2 2.4 2.6 

Fig. 4. Density evolution of the QCD susceptibilities (normalized to the vacuum value of the pseu- 
doscalar one) with = 8., m a — 850 MeV and C = 0.985. Dashed curve: pseudoscalar 
susceptibility. Full curve: Scalar susceptibility. Dotted curve: nuclear contribution to the scalar 
susceptibility. Dot-dashed curve: pion loop contribution to the scalar susceptibility. 



constrain the parameters of the ouj model. This, together with the information from 
spin-isospin physics, has allowed a successfull description of the binding properties 
of nuclear matter. 
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